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Abstract 

In a sample of 58 million J/ip events collected with the BES II detector, the process 
J/V' ^ JVc is observed in five different decay channels: 'yK'^K~7r'^Tr~ , ^-k^ it~ tt'^ it~ ^ 
^K^K^-K^ (with Kg — TT+TT"), (with ^ — K'^K~) and ^pp. Prom a combined 
fit of all five channels, we determine the mass and full- width of r]c to be m^^^ = 
2m. b ± 1.0 (stat.) ± 1.2 (syst.) MeV/c^ and P^^ = 17.0 it 3.7 (stat.) ± 7.4 (syst.) 
MeV/c^. 

PACS: 13.25.Gv, 14.40.Gx, 13.40.Hq 
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Since 1980, numerous efforts fiave been made to determine tlic mass and full 
width of the 1]^. [1-11]. However, from a theoretical point of view, the ac- 
curacies of these experimental measurements are still not sufficient. For in- 
stance, in order to calculate the strength of the spin-spin interaction term in 
non-relativistic potential models, it is necessary to know precisely the mass 
difference between the J/'0(1 ) and ric{0 ^). While the mass of the J/ijj is 
determined with high accuracy to be 3096.88±0.04 MeV/c^, the rjc mass is 
measured with much less accuracy to be 2979.7 ± 1.5 MeV/c^, an average 
by the Particle Data Group (PDG) [12] of 10 measurements with an internal 
confidence level of only 0.001. Different measurements of the full width of the 
Tjc also have poor internal consistency. The PDG [12] determines an average 
value for the r]c full width of 16.0l3;2 MeV/c^ from six experiments, whose 
experimental results vary from 7 MeV/c^ to 27 MeV/c^, with large errors. 
Such an accuracy is inadequate for some studies of charmonium physics [13] 
and additional, more precise measurements of both m^^ and P^^ are needed. 

The rjc mass and width have been measured previously by the BES collab- 
oration with data samples of 3.79 million ip{2S) events [9] and 7.8 million 
J/ijj events [11] collected with the BES I detector [14]. In the latter sam- 
ple, the process J/ip ^yrjc was observed in five different rjc decay chan- 
nels: K~^K~iT~^TT~ , 7r"'"7r~7r"'"7r~, K^Kgir^ (with Kg 7r"'"7r~), 00 (with 
K^K^) and K^K~ti^ , and the mass of the r/c was determined to be 2976.6 ± 
2.9 (stat.)±1.3 (syst.) McV/c^. Combined with the results from ■0(25') — > 7?7c, 
the mass and the full width of rjc were determined to be m^^ = 2976.3 ± 
2.3 (stat.) ± 1.2 (syst.) MeV/c^ and P^^ = 11.0 ± 8.1 (stat.) ± 4.1 (syst.) 
MeV/c^. 

In this paper we present results with much higher statistics using a recent 
sample of 58 million J/ip events obtained with the upgraded BESII detector 
[15]. The upgrade from BES I to BES II includes the replacement of the in- 
ner drift chamber with a straw-tube vertex chamber (VC), composed of 12 
tracking layers arranged around a beryllium beam pipe and with a spatial 
resolution of about 90 /xm; a new barrel time-of-fiight counter (BTOF) with a 
time resolution of 180 ps; and a new main drift chamber (MDC), which has 10 
tracking layers providing a dE/dx resolution of o^Ejdx — 8.4% and a momen- 
tum resolution of Up/p = 1.7%-\/l (p in GeV) for charged tracks. These 
upgrades augment the pre-existing calorimeter and muon tracking systems. 
The barrel shower counter (BSC), which covers 80% of 47r solid angle, has an 
energy resolution oi ge/ E = 22%/\/E (E in GeV) and a spatial resolution of 
7.9 mrad in and 2.3 cm in z. The /i identification system consists of three 
double layers of proportional tubes interspersed in the iron flux return of the 
magnet. They provide coordinate measurements along the muon trajectories 
with resolutions of 3 cm and 5.5 cm in and z, respectively. 

The 77c mass and width are measured using the reactions J/ip ^ ■jrjc', rjc 
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K+K-n+Ti-, 7r+7r-7r+7r-, K^K^n^ (with Tr+TT"), 00 (with ^ /sT+X-) 

and pp. Event selection criteria for each channel are described in detail in our 
previous papers [16,17,18]. Here we repeat only the essential information and 
emphasize those considerations that are unique to the m^^ and F^^ measure- 
ments. 

Candidate events are required to have the correct number of charged tracks 
for a given hypothesis. Each track must be well fit to a helix in the polar 
angle range | cos6'| < 0.84 and have a transverse momentum above 60 MeV/c. 
For the decay channels J/V' — > ^K'^ K~t:'^t^'' , J/'^ — > 77r''"7r~7r+7r~, J/V' — > 
^K'^T^'^T^^'K' and J/ip ^ "fpp, at least one photon with energy Ej > 30 MeV 
is required in the barrel shower counter 

Events are kinematically fitted with four constraints (4C) to the hypothe- 
ses: J/'^ — > "yK'^K'Tr'^Tr" , J/ip ^ 77r+7r~7r''"7r~, J/ip ^ 7X^71^71 '•"tt", and 
J/ijj ^ 7pp. A one-constraint (IC) fit is performed for the J/ip ^ '-fmissK'^ K~ K' 
hypothesis, where 'jmiss indicates that this photon is not detected. Events with 
a l^ss than 40.0 for a particular channel are selected. 

In order to remove backgrounds from non-radiative decay channels, all se- 
lected events arc subjected to a kinematic fit with four constraints to the hy- 
potheses: J/ip K~^K~7i~^'n'^ , J/ip — > 7r"'"7r~7r"'"7r~ and J/ip K^tt^tt^ti' . 
Backgrounds from the J/ijj peak are removed by requiring that x^{J/^ 
K+K-n+n-) > 20.0 (for K+K-n+n-); x^J/^ ^ n+n-n+n-) > 10.0 (for 
Tr+TT-TT+TT-) and x2(J/V' ^ K^n^n+n-) > 10.0 (for K^K'^s^^). For the 
J/ip ^ ^pp channel, we require that the opening angle of the two charged 
tracks is smaller than 179°. A detailed Monte Carlo simulation shows that 
these cuts, referred to below as the J/ip veto, do not distort the invariant 
mass distributions around the rjc signal peak. 

Two additional variables are used to reject events with wrong final state assign- 
ments. The first variable, |f/miss| = \Emiss — Pmiss\, is used to reject events with 
multi-photons and misidentified charged particles. Here, E^iss and Pmiss are, 
respectively, the missing energy and momentum calculated using measured 
quantities for charged tracks. A second variable, P^^ — 'i\Pjniss\'^ sm^{9tj/2), 
where Ofy is the angle between the missing momentum and the photon di- 
rection, is used to reduce backgrounds from vr'^'s. The specific values of the 
selection requirements for these two kinematic variables are summarized in 
Table 1. Additional requirements to remove backgrounds from a few specific 
channels are summarized in Table 2. 

For J/%1) ^ 'y7i~^n~n~^7i~ candidate events with more than one 7, we suppress 
7r° background by requiring that | M(7i72) — M{7r^) \> 60 MeV/c^ if Pmiss is 
in the same plane as the two photons 71 and 72, i.e. Pmiss ■ {i^^i x ''^72) < 0.15. 
Here, Pmiss is the unit vector of the missing momentum for all charged tracks; 
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Table 1 



Cuts imposed on |C/mjss| and for event selection. 



mode {J/ip —>■ ) 


1 Umiss 1 (GeV/c2) 


Pl^ [(GeV/c)2) 




< 0.15 


< 0.002 


77r"'"7r~7r+7r~ 


< 0.10 


< 0.0015 


7if±K§7rT (7if±7rT7r+7r-) 




< 0.003 




< 0.15 


< 0.003 



Table 2 

Cuts to remove backgrounds from specific channels. 



mode {J/ip iX ) 


cut 


background 




|M^+^-^o - M^l > 40 MeV/c^ 


J/ip LoK+K- 




\Mk+k- - > 20 MeV/c^ 


J/tp ^ (j)n^7r~ 


77r"'"7r~7r"'"7r~ 


|M^+^-^o - M^l > 40 MeV/c2 




77r"'"7r'~7r"'"7r~ 


|M^+^- - Mj^o \ > 25 MeV/c2 





r^i and r^2 arc unit vectors for the 71 and 72 directions determined from the 
BSC; and M(7i72) is the invariant mass of 7172- When we calculate M(7i72), 
it is assumed that the missing particle decays to 71 and 72, and M(7i72) can 

— * 

be obtained by using P„wss and the angles between Pmiss and the 7 direc- 
tion. The advantage of this technique is that it uses the momenta of charged 
tracks measured by the MDC, which has good momentum resolution, and is 
independent of the photon energy measurement. For J/ip ^ ■jK'^K~7r'^Tr~ , 
-fK^Kl-K^, and 7pp, we require that | M(7i72) - M(7r°) | > 50 MeV/c^ 
when Pmiss ■ (j^Ti X r:^2) < 0.14. 

For the K'^K^t:^ (with — >■ 7r"'"7r~) channel, the 7r"'"7r~ invariant mass for 
the K'g candidate is required to be within 25 MeV/c^ of the mass. For the 
00 (with (f) — > K^K^) channel, the invariant masses of both candidate 0's, 
corresponding to K^K~ pairs, are required to be within 20 MeV/c^ of the 
mass. 

After the event selection, the invariant mass spectra for the individual decay 
modes are obtained, as shown in Fig. 1. An unbinned maximum likelihood fit 
using MINUIT [19] is performed for all five channels simultaneously, with the 
fitting function for a given channel i given by 

fi{m) = ai[BW{M, F, m) ® GS{m, ai)] x EFFi{m) + (1 - ai)BGi{m), 

where M and F are the mass and width of the r]c, respectively, (Xj is the mass 
resolution in the rjc region, BW is a Breit-Wigner function describing the r]c 
signal, EFFi is an efficiency correction function, and BGi is a second-order 
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polynomial function describing the background shape. In order to include the 
experimental resolution, the BW function is folded with a Gaussian resolution 
function GS with the resolution Uj fixed at a value determined from the Monte 
Carlo simulation. The parameters M and F and the coefficients of the poly- 
nomial function, Cj, are determined from the fit. The log likelihood function 
for the channel i is given by 

jyevent 

i 

5"^ = -InLj = -ln( JJ fi{mj)), 
where Nf^^^^ is the total number of events. The overall log likelihood function, 

i=l 

is minimized to obtain the fitting results from the five channels simultaneously. 
The fit result is shown in Fig. 1, and the fitted ric mass and width are deter- 
mined to be m^c = 2977.5 ± 1.0 MeV/c^ and F^^ = 17.0 ± 3.7 MeV/c^. The 
background in Fig. 1(b), (d) and (e) can also be fitted with a linear polynomial 
function, with results that are almost the same. 

Systematic errors in determining the r]c mass and width originate mainly from 
the mass-scale calibration, background shape, fitting range, difference between 
data and Monte Carlo simulation, J/ip veto, and uncertainties associated with 
the mass resolution. We use 1.5 million ip{2S) data collected during the J/ip 
run to check the mass-scale calibration. The measured Xc2 mass is 3555.2 ± 1.4 
MeV/c^ from decays ^(25) ^ 'Jtt+tt-tt+tt- and 3560.2 ± 6.0 MeV/c^ from 
decays ■0(25') jK'^K~7r'^7r~, respectively. The combined weighted average 
is m^^2 = 3555.5 ± 1.4 MeV/c^, a difference of 0.7 ± 1.4 MeV/c^ from the 
world average obtained by the PDG [12]. In addition, we measured masses of 
the K^, (f) and A from the 58 million J/ip data sample to check the mass-scale 
calibration. Results of the masses and mass differences with PDG values [12] 
are given in Table 3. The systematic error on the overall mass scale is estimated 
to be 0.8 MeV/c^. 

Table 3 



Comparison of Kg, (p and A mass peak positions. 





i^o (MeV/c2) 


(MeV/c^) 


A (MeV/c2) 


our measurements 


496.9 ± 0.1 


1019.6 ±0.1 


1115.3 ±0.1 


PDG values 


497.67 ± 0.03 


1019.417 ± 0.014 


1115.683 ±0.006 


AM 


-0.8 ±0.1 


0.2 ±0.1 


-0.4 ±0.1 



Table 4 summarizes all contributions to the systematic error of the mass and 
full- width of the rjc- The effect of the background shape is studied by using 
a third-order polynomial function instead of a second-order one. The upper 
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Fig. 1. The invariant mass distributions in the r/c region for channels (a) 
mK+K-K+-K-: (b) "T-^+Tr-vr+Tr-, (c) mj^±f^o^^, (d) and (e) rupp. 



fitting bound is checked by changing it from 3.05 to 3.07 GeV/c^, and the 
J/ip veto is removed from the event selection. Contributions from differences 
of the detection efficiency between data and Monte Carlo simulation, as well 
as uncertainties of the detector mass resolution, are also listed in Table 4. 
Assuming no correlations among the above factors, the total systematic error 
on the mass and width are determined to be 1.2 MeV/c^ and 7.4 MeV/c^, 
respectively, by a quadratic sum of all contributions. 

In summary, we used 58 million J/ip events collected by the BES II detector 
to measure the mass and full width of the t]c in five different decay modes. 
They are determined to be m^^ = 2977.5 ± 1.0 (stat.) ± 1.2 (syst.) MeV/c^ 
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and r,,,^ = 17.0 ± 3.7 (stat.) ± 7.4 (syst.) MeV/c^. Figure 2 shows the BES 
results together with previously reported measurements. It can be seen that 
the r]c mass and width measurement from BES II are in good agreement with 
the PDG averages. 

Table 4 



Sources of systematic errors. 



sources 


error on mass 


error on width 




(MeV/c2) 


(MeV/c2) 


mass scale calibration 


0.8 




background shape 


0.0 


3.6 


fitting range 


0.3 


2.2 


detection efficiency difference: data vs. MC 


0.5 


0.3 


J/ip veto 


0.7 


5.6 


uncertainties of experimental mass resolution 


0.0 


2.4 


total systematic error 


1.2 


7.4 



PDG Value 2979.7± 1 .5MeV 

J 

^ CLE02 

^*^m DM2 

I • 1 ^ % MARK3(gO) 

»>-H SPEC 

^ MARK3(86) 

1 CBAL 

1 ^« 1 MARK2 

^ , ^ , , ^ 

2960 2970 2980 2990 3000 3010 MeV 




40 MeV 



Fig. 2. Mass and full-width measurements of the rjc meson. 
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